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Abstract. We present abundance measurements, based on high-resolution optical spectroscopic data obtained 
with the Ultraviolet and Visual Echelle Spectrograph mounted on Kueyen (Very Large Telescope UT2), for ten 
stars in the extended blue horizontal branch of the Galactic globular cluster NGC 1904 (M79). In agreement 
with previous findings for other clusters, we obtain normal abundances for stars cooler than T e g ~ 11000 K, 
and largely anomalous abundances for hotter stars: large He depletions, and overabundances of Fe, Ti, Cr, P and 
Mn. The abundances of Mg, Si and Ca are roughly normal, in the hot stars as well as in the cooler ones. This 
abundance pattern can be attributed to the onset of diffusion and to radiation pressure in the stable atmospheres 
of hot horizontal branch stars. 
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1. Introduction 

Sandage & Wildey (1967) and van den Bergh (1967) first 
noticed that the distribution of stars along the horizon- 
tal branch (HB) of globular clusters (GCs) cannot be de- 
scribed as a simple uniparametric function of metal abun- 
dance. In spite of intensive efforts since then, the addi- 
tional (second) parameter has not yet been clearly identi- 
fied; while numerous candidates have been proposed, none 
of them is sufficient to justify the whole pattern of obser- 
vations. The issue is likely complicated by several factors. 
In particular, the distribution of stars along the warmest 
side of the horizontal branch (i.e. stars bluer than the 
RR Lyrae instability strip, hence BHB) is very complex, 
exhibiting a large variety of shapes, from short, stubby 
horizontal branches (like e.g. in NGC 6397), to very long 
blue tails (like e.g. in NGC 1904). 

The distribution of stars along the HB of a GC is es- 
sentially determined by their masses: the smaller the mass, 
the bluer its colour. Since within a GC, stars have virtually 
the same age and chemical composition and hence possi- 
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bly the same mass at the turn-off (see however D'Antona 
et al. 2002), the colour of stars along the HB should be 
determined by the amount of mass lost by each star be- 
fore reaching the HB (the larger the mass lost, the bluer 
the colour on the HB). While it seems quite reasonable to 
assume that mass loss occurs more easily in extended and 
luminous red giants, the mechanisms driving mass loss are 
far from being well understood, in particular in the case of 
metal-poor, GC stars. Any possible observation relevant 
to this topic might help to constrain these mechanisms. 

Core rotation (Mengel & Gross 1976; Sweigart & 
Catelan 1998) is a reasonable candidate for enhanced mass 
loss along the red giant branch (RGB). In fact, core rota- 
tion reduces pressure in the helium cores, delaying the he- 
lium flash that terminates evolution along the RGB. Stars 
with rotating cores may then evolve to brighter luminosi- 
ties, and have the possibility to lose more mass. We would 
then expect that descendants of RGB stars with fast core 
rotation would inhabit the bluest part of the HB. Some 
support for this hypothesis was given in the '80s by ob- 
servation of rotating BHB stars in M13 (a cluster with 
a bluer than expected HB), but not in M3, by Peterson 
et al. (1983) and Peterson (1983, 1985a, 1985b). In recent 
years, with the advent of 8-10 m class telescopes, this issue 
has been more extensively revisited by Behr et al. (1999, 
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2000a, 2000b; revised in Behr 2003a, 2003b) and more 
recently by our group (Recio-Blanco et al. 2002). These 
additional observations revealed a much more complicated 
pattern than expected. Fast rotating stars are present also 
among field HB stars (i.e., this is not an effect limited to 
globular cluster stars, like the O-Na anticorrelation, see 
e.g. Kraft 1994), and are present in most clusters (even 
M3, according to Behr 2003a). They are however found 
only among the stars cooler than ~ 11 500 K. This tem- 
perature value plays an important role in the HB picture: 
stars warmer than this (and cooler than about 20 000 K) 
have anomalously bright u magnitudes in Stromgren pho- 
tometry (Grundahl et al. 1998, 1999); also, gravities de- 
rived for these stars from the wings of Balmer and he- 
lium lines are systematically lower than predicted by evo- 
lutionary models (Moehler et al. 1995, 2003). As shown 
by several authors (Grundahl et al. 1999; Moehler et al. 
1999; Hui-Bon-Hoa et al. 2000), this can be attributed 
to the onset of strong effects of levitation due to micro- 
scopic diffusion and radiation pressure that can be seen in 
stars lacking significant sub-atmospheric convection zones 
(Sweigart 2002). When microscopic diffusion and radia- 
tion pressure are not overwhelmed by convection, the at- 
mospheres of warm stars are expected to be strongly de- 
pleted in He (Greenstein et al. 1967; Michaud et al. 1983), 
and to exhibit large overabundances of elements like Fc, 
Ti and P (Turcotte et al. 1998; Richer et al. 2000). All 
these features are indeed observed (Glaspey et al. 1989; 
Moehler et al. 1999; Behr 2003a). 

It should be noted that most of these results - in partic- 
ular those related to the chemical composition - have been 
obtained only for a few GCs, mainly by a single group. 
It would be useful to extend such observations to more 
clusters, performing an independent analysis. In this pa- 
per we present a new analysis of the chemical composi- 
tion of 10 BHB stars in the southern cluster NGC 1904 
(M79). NGC 1904 is a cluster of intermediate metal abun- 
dance ([Fe/H]=-1.59 according to Kraft & Ivans 2003, us- 
ing Kurucz 1994 models), with an extended very blue HB: 
these properties are similar to those of M13, one of the 
clusters studied by Behr (2003b). The observational ma- 
terial obtained by us for NGC 1904 was originally intended 
only for measurements of the rotational velocities: hence 
the spectra have a low signal-to-noise ratio (S/N), al- 
though high resolution. However, we will show that abun- 
dances good enough for the purposes of the present dis- 
cussion could be obtained even from this material. 

2. Observations and reduction 

A total of 21 stars were observed in NGC 1904; however, 
spectra for only 10 of these stars were deemed of good 
enough quality for abundance analysis, the criteria being 
the S/N of the spectra and the rotational velocity of the 
stars (only slow rotators were considered) . We selected our 
targets from the F439W and F555W HST-WFPC2 pho- 
tometry by Piotto et al. (2002). In addition, we have also 
identified the NGC 1904 targets in the Stromgren u, y pho- 



tometry by Grundahl et al. (1999) and the Johnson U, V 
photometry by Momany et al. (2004). In general, targets 
lie in the low-density outskirts of the GC, to avoid contam- 
ination from other stars. The spectra were obtained dur- 
ing 2 observing runs: July 30- August 2 2000 and January 
19-23 2001. 

The observations were carried out with the UVES 
echelle spectrograph, mounted on the Unit 2 (Kueyen) of 
the VLT, and the 2K x 4K, 15 ^m pixel size blue CCD, 
with a readout noise of 3.90 e~ and a conversion factor of 
2.04 e"/ADU. The UVES blue arm, with a spectral cov- 

o 

erage in the 3730-5000 A range, was used. Combined with 
a slit width of 1.0 arcsec, we achieved a resolution of R~ 

40 000 (SX ~ 0.1 A, Sv ~ 7.5 km/s, see the UVES user 
manual by Kaufer et al. 2003). 

The exposure times ranged from 800 s (for V ~ 16) to 
4000 s (for V ~ 17.5). We generally limited individual ex- 
posure times to 1500 or 2000 seconds, to minimize cosmic 
ray accumulation, and then coadded individual frames ob- 
tained for each star. The typical signal to noise ratios of 
the coadded spectra are about 20 to 30 per resolution ele- 
ment, but in a few cases they are as good as S/N ^40-60. 
See Fig. for a visual example of the quality of the target 
spectra. 

The spectra were extracted using standard IRAF pro- 
cedures. For the wavelength calibration we fitted 3 rd order 
polynomials to the dispersion relations of the ThAr cali- 

_ . o 

bration spectra, which resulted in residuals of < 3T0 A- 
Finally, each spectrum was normalized using a 5th-degree 
polynomial fitting. 



3. Abundance measurements 

We have performed an abundance analysis of 10 stars in 
M79. The abundance measurements were performed us- 
ing the program WIDTH3, developed by R. G. Gratton 
and adapted by D. Fabbian to temperatures up to T e g ~ 
20 000 K (D. Fabbian graduate thesis). The procedure es- 
tablishes the abundances of chemical species by repro- 
ducing the observed equivalent widths. Once a starting 
set of values for the effective temperature, surface grav- 
ity and model metal abundance is derived, an appro- 
priate model for each star is obtained from the grid of 
model atmospheres by Kurucz (1994) by interpolating 
linearly in temperature and logarithmically to get the 
Rosseland opacity, electronic and gaseous pressure, and 
density. Continuum opacity was obtained taking into ac- 
count all important continuum opacity sources for stars 
as hot as T e g — 20 000 K. Collisional damping constants 
were computed with the Unsold formula. The equation of 
transfer was then integrated through the atmosphere at 
different wavelengths along the line profile and theoreti- 
cal equivalent widths of the lines computed and compared 
with the observed ones. 
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Fig. 1. Spectral sections of the stars 392 (top panel), and 209 (bottom panel) 



3.1. Oscillator strengths, line selection and solar 
abundances 

Only lines free of blends were considered in the analysis. 
The line list was extracted from the solar spectrum tables 
by Moore et al. (1966), and the atomic line lists by Hambly 
et al. (1997) and Kurucz & Bell (1995), more appropriate 
for warmer stars. 

For the determination of the oscillator strengths (gjs), 
laboratory values were considered whenever possible. For 
Fe I lines they were taken from papers of the Oxford group 
(Simmons & Blackwell 1982) and Bard & Kock (1994). 
The gjs for the Fell lines were taken from Heise & Kock 
(1990) and Biemont et al. (1991), while gjs for the Mgl 
lines came from Gratton (2003). For the rest of the lines, 
the gjs were those provided by Kurucz (1994) and the 
NIST atomic spectra database. 

The solar abundances were taken from Grevesse & 
Sauval (1998). 



3.2. Equivalent Widths 

Equivalent Widths (EWs) were measured on the unidi- 
mcnsional, extracted spectral orders using an automatic 
routine within the ISA package, prepared by R. Gratton 
(see Bragaglia et al. 2001). The routine works as follows. 
First, a local continuum level is determined for each line 
by an iterative clipping average of the 200 spectral points 
centered on the line to be measured. Second, the equiva- 
lent width for each line from an extensive list is tentatively 
measured using a Gaussian fitting routine. Measures are 
rejected according to several criteria (e.g., if the central 
wavelength does not agree with that expected from a pre- 
liminary measure of the geocentric radial velocity or if the 
lines are too broad or too narrow). Third, the measured 
lines are used to draw a relationship between equivalent 
width and FWHM. Fourth, this relation is used to mea- 
sure equivalent widths again, using a second Gaussian fit- 
ting routine that has only one free parameter (the equiv- 
alent width), since the central wavelength is fixed by the 
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Table 1. The measured equivalent widths and corre- 
sponding uncertainties used in our abundance analysis. 
The table is available only in electronic form at the CDS. 
Column 1 lists the chemical elements of interest with their 
ionization stage; Column 2 gives the corresponding wave- 
lengths (in A); excitation potentials (in eV) and log gjs are 
given in Column 3 and 4 respectively; the measured equiv- 
alent widths and corresponding uncertainties (in m A) for 
the 10 stars are listed in the next columns, with a few 
values being upper limits 



Table 2. Abundance errors for the star 392, due to un- 
certainties a(AEW) in equivalent width, a (AT) in effec- 
tive temperature, a(Ag) in logg, c(A£) in £ and cr(Az)in 
[A/H] 




Element 


a(AEW) 


a(AT) 


a(Ag) 


a(A0 


a(Az) 




±10.4 


±200 


±0.1 


±1 


±0.1 




(ml) 


(K) 


dex 


km/s 


dex 


He i 


0.3 


0.11 


0.06 


0.02 


0.09 


Mgn 


0.2 


0.02 


0.01 


0.08 


0.00 


Sin 


0.2 


0.02 


0.03 


0.03 


0.01 


Pn 


0.2 


0.02 


0.03 


0.05 


0.01 


Can 


0.4 


0.05 


0.04 


0.22 


0.01 


Tin 


0.2 


0.09 


0.00 


0.05 


0.02 


Cm 


0.2 


0.04 


0.02 


0.03 


0.02 


Mnii 


0.2 


0.03 


0.02 


0.02 


0.01 


Fei 


0.2 


0.10 


0.03 


0.05 


0.04 


Fen 


0.2 


0.02 


0.03 


0.09 


0.01 



measured geocentric radial velocity and the FWHM is 
fixed by the relationship between equivalent width and 
FWHM. Again measures are rejected if residuals are too 
large, asymmetric, etc. 

But for a few exceptions (Mg n and Ca II lines having 
EW >100 mA ), only lines with 10< EW <100 mA were 
finally considered. Lines with EW <30 mA were only used 
for the star 392, the one with the best spectrum (see sec- 
tion 4.2 and Figure^). The number of measured lines de- 
pends on the S/N and on the star's metallicity. Generally, 
around 15 Fei lines and 20 Fen were measured. Measured 
EWs for individual stars are listed in Tabled on ly avail- 
able in electronic form. 

Due to the quality of the spectra, the EW uncertanties 
were the most important source of errors in the abundance 
determinations. This error was evaluated by comparing 
the equivalent width measurements of the stars 209 and 
281. These stars have almost equal atmospheric param- 
eters, with differences just slightly larger than the inter- 
nal errors. In addition, they do not present signs of diffu- 
sive effects in their atmospheres. An analysis of 22 lines 
measured in both spectra indicates that the mean differ- 
ence between the two equivalent width measurements is 
— 1.1 mA (in the sense 281-209), with a standard deviation 
of 14.7 mA. Therefore, assuming that equivalent width er- 
rors for both stars are the same, the error for each mea- 
surement is 14.7/v / 2=10.4 mA. If we compare this result 



Fig. 2. Colour magnitude diagram for the target stars. 
Filled symbols are stars for which abundance analysis was 
done; open symbols are those stars for which only rota- 
tional velocities were obtained. Left panel: U — (U — V) 
colour magnitude diagram from Johnson photometry by 
Momany et al. (2004). Right panel: u,(u — y) diagram from 
Stromgren photometry by Grundahl et al. (1999). ZAHB 
models by Cassisi et al. (1999) are plotted as continuous 
lines. The location of the luminosity jump is marked with 
an arrow in both CMDs 

with the value calculated through the Cayrel (1988) for- 
mula: 

A(EW) = l.6 v / (wdx)/(S/N) 

with w being the full width at half maximun (FWHM) 
typical of the lines, in this case ^0.12 A, and dx being 
the pixel size, about 0.03 A, we get A(EW) ~ 10 mA 
(assuming a S/N ~10 per pixel), which is more or less 
what we obtained comparing the measurements for the 
stars 209 and 281. This result implies that identification 
of the correct continuum level (that is neglected in the 
Cayrel formula) is not an issue here, as expected since the 
spectra of BHB stars have very few absorption lines. 

Consequently, we calculated the abundance errors, de- 
rived from an equivalent width uncertanty of 10.4 mA. 
The resulting errors for the star 392, a(AEW), are shown 
in Column 2 of Tabled 

3.3. Atmospheric parameters 

Model atmospheres appropriate for each star were ex- 
tracted from the grid of Kurucz (1994). Figure [21 shows 
the projections of the analyzed HB stars into the corre- 
sponding zero age horizontal branches (ZAHBs) for the 
two adopted photometries. The position of the Grundahl 
et al. (1999) luminosity jump is marked with a vertical 
arrow. This feature clearly affects the fit of the models 
to the observed HBs, especially in the u,(u — y) CMD. 
T e ff values were determined from the position of the stars 
in the u,(u - y) CMD by Grundahl et al. (1999), and 
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Table 3. Photometry for all the target stars in NGC 1904. The temperature values derived from the u,(u — y) Grundahl 
et al. (1999) CMD, T e s g , the corresponding corrected effective temperature values T e g gc and the temperature derived 
from the Momany et al. (2004) CMD, T eSm 
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17.55 
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the position in a Johnson U,(U — V) CMD by Momany 
et al. (2004), using the ZAHB models by Cassisi et al. 
(1999), see Fig- EI Those models were transformed from 
the theoretical plane to the observational one using the 
colour-effective temperature relation provided by Castelli 
et al. (1997a, b). Hereafter, we adopt a reddening value 
of E(B - V) = 0.01, from the Harris (1996) compila- 
tion. The value that could be obtained using the maps 
of Schlegel et al. (1998) based on COBE-DIRBE data is 
E(B — V) — 0.031. If we had rather adopted this last value 
of the interstellar reddening, the temperatures would have 
been ~ 200 K warmer; this is less than the internal errors. 

Table |3 presents the photometric parameters for all 21 
stars in NGC 1904 for which we obtained spectra, with 
the corresponding temperature values derived from the 
Grundahl et al. (1999) u,(u - y) CMD and the Momany 
et al. (2004) CMD, respectively labelled T cffg and T cffm . 
Only 10 of these stars, selected among the slow rotating 
ones, were used in the abundance analysis. 

To derive the final effective temperature of the pro- 
gramme stars, we calculate the mean of Momany's tem- 
peratures and of Grundahl's temperatures after transpos- 
ing these last to Momany's scale. The final error is the 
scatter of the temperatures involved in the mean. We de- 
cided to convert the T e g E values to Momany's scale mainly 
because the T c g g scale presented a number of problems. 
For instance, the T e g g temperature for the coolest stars 
would place them inside the RR Lyrae instability strip, 
and the corresponding metallicity for such a cool tem- 
perature would be lower than what is reported for the 
RGB stars in the cluster. On the other hand, apart from 



Table 4. Atmospheric parameters for stars in NGC 1904 



star 




logs 


[A/H] 






(K) 






(km/s) 


209 


8469±200 


3.2±0.2 


-1.4 


2.4±1.0 


281 


8925±200 


3.3±0.2 


-1.4 


2.3±1.0 


354 


10409±200 


3.7±0.1 


-1.4 


2.0±1.0 


389 


11148±200 


3.8±0.1 


0.0 


O.Oil.O 


363 


11681±300 


3.9±0.1 


0.0 


O.Oil.O 


392 


11700±200 


3.9±0.1 


0.0 


O.Oil.O 


434 


11823±200 


3.9±0.1 


0.0 


O.Oil.O 


469 


12123±200 


4.0±0.1 


0.0 


O.Oil.O 


535 


13070±300 


4.1±0.1 


0.0 


O.Oil.O 


489 


13327±200 


4.2±0.1 


0.0 


O.Oil.O 



this shift to cooler temperatures, probably due to the 
ZAHB model fit, all of our targets were well identified 
in the Grundahl et al. photometry, while 3 of the stars 
in Momany's photometry were blended. For this reason, 
and to be sure of an internal consistency in all the tem- 
perature measurements, we decided to transform the T ff g 
values into Momany's scale, and then calculate the mean, 
for each star, of the new T c ff g obtained and the T g m value. 

The relation between T e g g and T e ff m was found by av- 
eraging the best fit relations obtained assuming the two 
quantities as independent variables, because internal er- 
rors in T e ff m and T e ff g were expected to be comparable. 
The final relation is: 

Tefigc = (0.731 ± 0.039)T cffg -I- (3055 ± 357) K (1) 
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Fig. 3. Runs of T e ff against the dereddened Johnson B — 
V color (left panel) and absolute My magnitude (right 
panel) for our program stars in NGC 1904 (filled symbols), 
and for the stars in M13 analyzed by Behr (2003a) (open 
symbols) 



where T c ffgc is the corrected effective temperature de- 
rived from the Stromgren u, y photometry by Grundahl 
et al. (1999). Internal errors in these temperatures may be 
obtained by comparing the differences between the tem- 
peratures obtained from the two photometries. The r.m.s. 
average of these differences (taking into account the two 
degrees of freedom lost because we corrected the temper- 
atures obtained from GrundahPs photometry to those ob- 
tained from Momany's photometry) is 355 K; assuming 
that individual errors are similar in the two determina- 
tions, the expected internal errors are ±177 K for those 
stars with both determinations, and ±251 K for stars 
having only Grundahl's photometry. Hereafter, we will 
round these values to ±200 and ±300 K respectively 1 , 
see Table H 

In order to compare the effective temperatures adopted 
throughout this paper with those adopted in other abun- 
dance analyses of BHB stars, we plotted (see Fig. [3] left 
panel) the runs of T e g against the dereddened Johnson 
B ~ V color and of the absolute My magnitude for both 
our program stars in NGC 1904 and the stars in M13 
analyzed by Behr (2003a). The reddenings adopted here 
are those from the Harris (1996) compilation. There is a 
reasonable agreement between the two sets for the coolest 
stars, while our temperatures for the hottest ones are lower 
than those obtained by Behr for stars of similar B — V 
colours. This difference is likely due to some differences 
in the colour calibration; in fact, the agreement between 
the two sets is excellent if absolute My magnitudes are 
considered (see Fig.0 right panel). 



1 Systematic errors have of course larger values than these 
internal errors 



Fig. 4. Comparison between the gravities adopted in this 
paper (log pet) and those obtained from the location of the 
stars in the colour-magnitude diagram, using theoretical 
models by Cassisi et al. (1999: loggtco) 

Gravities are not well costrained by our spectra: in fact 
equilibrium of ionization is subject to possible departures 
from LTE, and the Balmer lines are too broad for reliable 
determinations of their profile from our Echelle spectra. 
We have then derived a mean relation between logT ff 
and gravity (logg) from Behr et al. (1999) measurements 
of blue HB stars in M13 that were expected to be very 
similar to those in NGC 1904 on the basis of the colour- 
magnitude diagram: 

log<? = 4.83 -log T eff -15.74 (2) 

Figure 0] compares the gravities obtained using this 
method (and used by us in the abundance determina- 
tions) , with those that could be derived from photometric 
angular diameter. In the latter case, a distance modulus 
of (m-M)\/ = 15.59 (Harris, 1996, in the revised version 
of 2003) has been adopted for the cluster, and a reason- 
able mass for the stars was estimated from HB models by 
Cassisi et al. (1999) for the Johnson photometry, again 
using the location of the stars in the colour-magnitude di- 
agram. The agreement between these two values is very 
good, supporting the gravity estimates adopted through- 
out this paper. Errors in the gravities due to errors in the 
temperatures are of about 0.05 dex. 

Microturbulent velocities £ might be derived by elim- 
inating trends of the derived abundances with expected 
line strength for some given species. However, given the 
quality of our spectra, generally too few lines could be 
measured, and with too large scatter, to significantly con- 
strain microturbulent velocity. For the cooler stars, we re- 
lied on Behr's analysis of M13 stars. From his analysis, 
we derived the following relation between logT c g and the 
microturbulence velocity: 

i= -4.7 -log T eff + 20.9 km/s (3) 
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Uncertainties in these values for the microturbulcnt ve- 
locities can be obtained from the scatter of individual val- 
ues around this mean relation: this is ~ 1 km/s. For the 
warmer stars, we adopted no microturbulent velocity, as 
given by the analysis of star 392, which has the best spec- 
trum (see Section 4.2). 

Metallicities were obtained varying the metal abun- 
dance [A/H] of the model until it was close to the de- 
rived [Fe/H] value. The adopted atmospheric parameters 
are listed in Table The errors. 

Once T e fi, logg, £ and atmospheric [Fe/H] were deter- 
mined, we calculated the mean abundance and dispersion 
for each element. Since the abundances for each element 
depend upon the adopted parameters, we recalculated the 
abundances varying each of these parameters in turn, to 
test how much the abundances change. The resulting er- 
rors for the star 392 are shown in Tabled Columns 3 to 
6. 

4. Results 

4.1. Abundance anomalies 

In this section we present the results of the abundance 
analysis. Tables [S] and [fj] contain the abundance values de- 
rived for each of our target stars, both on an absolute 
scale, and relative to the solar abundances of Grevesse & 
Sauval (1998). The scatter of the measurements for the 
corresponding lines gives the quantity ojv, which is added 
quadratically to the error contributions in T c ft, g, £ and 
[A/H] described in the previous chapter, to get the final 
error for each element in each target star. 

4.2. Star 392 

We will first consider the analysis of star 392, whose spec- 
trum has the highest S/N (~ 60 per resolution element 
at a wavelength of 4240 A). For this star we were able to 
measure over 100 lines. 

Figure [S] shows the curve of growth of Fen. Clearly, 
many unsaturated lines were accurately measured. As 
shown in Figure no trend was detected in the num- 
ber of Fe atoms, log n(Fe), with the excitation potential, 
as expected for the excitation equilibrium. This test is 
very critical, given the large spread of excitation poten- 
tials (nearly 10 eV): even a modest error of 500 K in the 
effective temperature would result in a very significant 
trend in this plot. 

The large number of measured lines for Fen allowed 
an accurate estimate of the microturbulent velocity: best 
results were obtained setting it to zero (i.e., leaving only 
thermal broadening as an unsaturating factor). This result 
is not surprising in view of the fact that this star is so hot 
that the atmosphere should be in radiative equilibrium. 

With our value of the effective temperature and sur- 
face gravity, we got an excellent equilibrium of ioniza- 
tion for Fe: the mean abundance given by neutral lines 
is log n(Fe)=7.86 ± 0.03 (31 lines, r.m.s. of results for in- 




-2 -1 

log(A)+log(gf)-0*EP 



Fig. 5. Curve of growth for Fell lines in the spectrum of 
star 392 




log(X) + log(gf)-0*EP 

Fig. 6. Upper panel: run of the extracted abundances for 
individual Fen lines from the spectrum of star 392 with 
excitation potential (in EV). Lower panel: the same, but 
against expected line strength 



dividual lines of 0.18 dex), while that obtained from singly 
ionized lines is log n(Fe)=7.93 ± 0.03 (48 lines, r.m.s. of 
results for individual line of 0.20 dex). The offset between 
the abundances is only marginally larger than the sum of 
the internal errors; again, a small change of the effective 
temperatures and gravities (within the internal errors of 
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Table 5. Abundances for stars in NGC 1904. The superscripts next to the abundance determinations indicate the 
number of lines used in each case 



star 


log He 


log Mg 


log Si 


log P 


log Ca 


log Ti 


log Cr 


log Mn 


log Fe 


209 




6.4±0.3 2 


6.2±0.2 X 




s^io.s 1 


3.7±0.2 13 


4.4±0.2 1 


< 4.7 1 


6.03±0.20 y 


281 




6.7±0.3 2 


6.5±0.2 3 






3.9±0.2 12 


4.8±0.3 1 




6.14±0.20 7 


354 


ll.SiO.3 1 


6.2±0.3 2 


6.4±0.2 2 






4.3±0.2 5 






6.34±0.20 6 


389 




6.5±0.3 2 


7.3±0.3 1 


7.4±0.3 1 


s.eio.s 1 


4.9±0.2 8 






7.90±0.20 14 


363 


< 9.4 1 


6.2±0.3 2 






e.sio.s 1 


5.3±0.2 5 






8.00±0.20 13 


392 


lO.giO.3 1 


6.3±0.3 2 


5.2±0.2 2 


6.8±0.3 6 


5.7+0.3 1 


4.8±0.2 12 


4.9±0.3 3 


e^io^ 1 


7.93±0.20 48 


434 


< 9.3 1 


6.5±0.3 2 


6.1±0.3 2 




e.oio.3 1 


5.1±0.2 6 


6.5ZL0.5 1 




7.90±0.20 24 


469 


g.sio.3 1 


6.7±0.3 2 


e.eio^ 1 




6.0±0.2 2 


5.4±0.2 7 


5.8±0.5 2 


6.5+0.3 1 


8.11±0.20 28 


535 


g.sio.3 1 


6.2±0.3 2 


7.2±0.2 4 




5.8+0.3 1 


4.6±0.2 1 






7.90±0.20 3 


489 


g.sio.s 1 


6.6±0.3 2 


7.1±0.2 3 


6^0.3* 




5.6±0.2 2 


7.4±0.3 1 




7.87±0.20 21 


Table 6. Abundance relative to solar for stars in 


NGC 1904 










star 


[He/H] 


[Mg/H] 


[Si/H] 


[P/H] 


[Ca/H] 


[Ti/H] 


[Cr/H] 


[Mn/H] 


[Fe/H] 


209 




-1.1±0.2 


-1.3±0.2 




-0.6±0.3 


-1.3±0.2 


-1.3±0.2 


< -0.7 


-1.48±0.20 


281 




-0.8±0.2 


-1.0±0.2 






-1.1±0.2 


-0.9±0.3 




-1.37±0.20 


354 


+0.3±0.3 


-1.3±0.2 


-1.1±0.2 




-0.9±0.3 


-0.7±0.2 






-1.17±0.20 


389 




-1.0±0.2 


-0.2±0.3 


+1.9±0.3 


-0.7±0.3 


-0.1±0.2 






+0.39±0.20 


363 


< -1.6 


-1.3±0.2 






0.0±0.3 


+0.3±0.2 






+0.49±0.20 


392 


-0.1±0.3 


-1.2±0.2 


-2.3±0.2 


+1.3±0.3 


-0.6±0.3 


-0.2±0.2 


-0.8±0.3 


+0.8±0.2 


+0.42±0.20 


434 


< -1.7 


-1.0±0.2 


-1.4±0.2 




-0.3±0.3 


+0.1±0.2 


+0.8±0.5 




+0.39±0.20 


469 


-1.2±0.3 


-0.8±0.2 


-0.9±0.3 




-0.3±0.2 


+0.4±0.2 


+0.1±0.5 


+1.1±0.3 


+0.60±0.20 


535 


-1.2±0.3 


-1.3±0.2 


-0.3±0.2 




-0.5±0.3 


-0.4±0.2 






+0.39±0.20 


489 


-1.2±0.3 


-0.9±0.2 


-0.4±0.2 


+1.2±0.3 


-0.1±0.3 


+0.6±0.2 


+1.7±0.3 




+0.36±0.20 



both these quantities) would bring the two abundances in 
agreement with each other. 

This very good agreement supports our LTE analysis. 
Actually, such a good agreement was not foreseen, since 
we expected some overionization of Fe, leading to lower 
abundances from neutral lines in an LTE analysis. This 
result should be compared with appropriate predictions 
from statistical equilibrium calculations for these stars; 
however, this is beyond the scopes of the present analysis. 

4.3. Other stars 

The abundance values for the various elements for which 
we have data are plotted in Figures |H1 ED HOI EH In 
the upper panels, the abundances for individual stars in 
NGC 1904 are plotted as a function of T e g as log offsets 
from the solar abundances. The abundance values from 
Behr (2003a) for M13 are plotted for comparison in the 
lower panels. The horizontal lines represent the expected 
value of [metal/H] from solar for NGC 1904 and M13 in 
the Carretta & Gratton (1997) scale: at -1.37 and -1.39, 
respectively. For helium, the horizontal line represents the 
solar ratio. We used the abundances derived from lines of 
the dominant ionization stage of each element, to mini- 
mize the possibility of non-LTE effects. 

4.3.1. Iron, titanium and chromium 

Statistics of abundances for Fe are given in Table 
Abundances are obtained from both neutral and singly 



Table 7. Statistics of Fe Abundances 



star 




Fel 






Fen 






n 


<> 


(7 


n 


<> 


O" 


209 


29 


6.06 


0.19 


9 


6.03 


0.25 


281 


19 


6.47 


0.27 


7 


6.14 


0.26 


354 


2 


6.90 


0.18 


6 


6.34 


0.32 


389 


11 


7.66 


0.27 


14 


7.90 


0.31 


363 


8 


8.24 


0.47 


13 


8.00 


0.38 


392 


31 


7.86 


0.18 


48 


7.93 


0.20 


434 


9 


7.80 


0.44 


24 


7.90 


0.48 


469 


9 


7.87 


0.26 


28 


8.11 


0.34 


535 








3 


7.90 


0.06 


489 








21 


7.87 


0.46 



ionized lines, in a few cases from a total of more than 
30 lines. Figure displays the run of the differences 
between the average abundances obtained from neutral 
and singly ionized Fe lines as a function of temperature. 
Abundances from Fe I lines agree well with those from Fe n 
lines. However, given the possibility that small non-LTE 
effects are present, in the following discussion we consid- 
ered the abundances provided by singly ionized Fe, which 
is the dominant species and should then be less affected 
by model uncertainties. 

As can be seen from Figure |H1 the iron abundances 
for the three coolest stars are very close to the cluster 
metallicity from analysis of red giant branch stars. We 
obtain [Fe/H]=— 1.48, —1.37, and —1.17, yielding an av- 
erage value of [Fc/H]=— 1.34 ± 0.09. This value may be 
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8000 10* 1.2x10* 1.4x10* 

T eff 

Fig. 7. Differences between abundances for Fe I and Fe II 
lines as a function of temperature 



compared with those obtained from analysis of red giants: 
[Fe/H]=-1.37 on the Carretta & Gratton (1997) scale and 
the slightly lower value of [Fe/H]=-1.59 given by Kraft & 
Ivans (2003). 

On the other hand we find remarkable enhancements 
of iron and other metal species among the blue HB stars 
hotter than ~ 11 000 K. The coolest star showing signs of 
radiative levitation is star 389 at T cS = 11 148 ± 200 K, 
with a supra-solar value of [Fc/H]= +0.39 ± 0.20. Stars 
hotter than ^11000 K show a very similar iron content: 
the average value is [Fe/H] =+0.43 ± 0.05. The very small 
scatter (0.08 dex) may be justified by the internal errors in 
the atmospheric parameters. Both the threshold temper- 
ature at which the anomalous abundances appear, as well 
as the amplitude of the enhancement (about 1.8 dex) are 
similar to those found by Behr et al. (2003a) for various 
clusters. 

In a similar way, titanium is found a few tenths of a dex 
above the cluster baseline in the cooler stars ([Ti/H]=-1.3, 
-1.1, and -0.7). The average value of [Ti/Fe]=+0.31 agrees 
very well with the overabundance of this element as gen- 
erally found in cooler metal-poor stars. Titanium abun- 
dances then rise by a factor quite similar to that found for 
Fe in the hotter population. 

A similar trend is also obtained for chromium, again in 
good agreement with the results by Behr (2003a); however, 
Cr abundances are based on only two Cm lines. 
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Fig. 8. Abundance values for iron, titanium and 
chromium as a function of T e g as log offsets from the so- 
lar abundances. The abundance values from Behr et al. 
(2003a) for M13 are also plotted for comparison in the 
lower panels. The horizontal solid lines represent the ex- 
pected value of [metal/H] from solar for NGC 1904 and 
M13, in the Carretta & Gratton (1997) scale. The ar- 
rows at 11 000 K indicate the approximate position of 
Grundahl's photometric jump 

4.3.2. Helium 

The He abundances appear below the expected 2 solar 
He/H ratio in the stars with T D g > 11 000 K. However, at 
variance with the case of Fe, we found a large and likely 
significant star-to-star scatter in the He depletions, see 
Fig. The case of star 392 is very significant: the rela- 
tively high S/N allows clear detection of the He lines. The 
He abundance we derive from the 4471 A line is consistent 
with the solar value. On the other hand, no He line at all 
was detected in other stars having a similar temperature. 
Weak lines were detected in the spectra of the warmer 
stars, suggesting a large He depletion of at least an order 
of magnitude. 

4.3.3. Phosphorus and manganese 

Some metals, such as phosphorus and manganese (see 
Fig. I10p . display significantly larger enhancements than 
iron. We do not observe any Pn line among the cooler 
stars, but if we assume an appropriately-scaled solar com- 
position for these stars, then the supra-solar value of [P/H] 
~ +1.5 that we find for the hot stars implies an enhance- 

2 We neglect here the difference between the solar and cos- 
mological He abundance, much smaller than the observational 
errors of our He abundance determinations. 
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Fig. 11. Same as Figure|Hl but for magnesium, silicon and 
calcium 
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Fig. 10. Same as Figure HO but for phosphorus and man- 
ganese 

ment of ~3 orders of magnitude with respect to the clus- 
ter's adopted metallicity. It is worth noting that several 
P II lines were found in the field HB stars Fcigc 86 and 3 
Cen A by Sargent & Searle (1967) and Bidelman (1960), 
respectively, suggesting that the same mechanism may be 
at work in both cluster stars and field stars. On the other 
hand, Behr (2003a) found very similar P enhancements for 
hot HB stars in M13. Similar very large overabundances 
are also observed for manganese. 



4.3.4. Magnesium, silicon and calcium 

Figure ^2 displays the run of the abundances of magne- 
sium, silicon and calcium with temperature. The abun- 
dance of magnesium was obtained from the strong multi- 
plet of Mgn at 4481 A, that is partially resolved in our 
spectra; Ca abundance was derived from the resonance 
K-line. 

The abundances for magnesium, silicon and calcium, 
both below and above the critical temperature of 11 000 K, 
are consistent with very little or no enhancement, but with 
a large, likely real scatter for Si. This result agrees with 
that found by Behr (2003a) and with theoretical expecta- 
tions. 

5. Discussion and conclusions 

The abundance anomalies in NGC 1904 are likely to be 
due to the same diffusion processes that were invoked by 
Behr et al. (1999): radiative levitation of metals and grav- 
itational settling of helium, in the stable non-convective 
atmospheres of the hotter, higher gravity stars, as hypoth- 
esized by Michaud et al. (1983). In addition, the jump 
in the u,(u — y) diagram by Grundahl et al. (1999) for 
NGC 1904 is reported at logT eff = 4.06 ± 0.01, in very 
good agreement with the onset of the abundance anoma- 
lies observed here. 

An underabundance of helium on the BHB was ob- 
served in several previous instances (Baschek 1975; Heber 
1987; Glaspey et al. 1989; among others) and in fact ap- 
pears to be typical of stars of this type. Michaud et al. 
(1983), building on the original suggestion by Greenstein 
et al. (1967), explained these underabundances as a re- 
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suit of the gravitational settling of helium, which can take 
place if the outer atmosphere of the star is sufficiently sta- 
ble. Our current results confirm a distinct trend in [Hc/H] 
with T c g and logg along the HB. Helium depletion is ac- 
companied by photospheric enhancement of metals, since 
the same stable atmosphere that permits gravitational 
settling also allows the levitation of elements with large 
radiative cross sections. Elements which have sufficiently 
large cross-sections to the outgoing radiation field will ex- 
perience radiative accelerations greater than gravity and 
will diffuse upwards, enriching the photosphere. The gross 
behaviour of the diffusion models can be understood by 
considering the migration of chemical species in a stellar 
atmosphere (Gonzalez et al. 1995; Michaud et al. 1983 
and Seaton 1997). The relative amount of levitation ex- 
perienced by each element can be estimated by summing 
the predicted equivalent widths for each line over the vis- 
ible and near UV spectral range. The magnitudes of these 
net equivalent widths suggest that helium will experience 
only a weak levitation force, insufficient to oppose gravity, 
magnesium will be levitated by an intermediate amount, 
perhaps balanced by gravity, while iron will be strongly 
levitated. This is admittedly a crude approach to what we 
see in NGC 1904, and Behr et al. measurements in various 
clusters. 

The identification of the detected spectral lines was 
successful in most cases. A few weak lines were not iden- 
tified, which could be due to the low S/N of the spectra. 
Since the abundances of metals we derived are very pe- 
culiar, it is possible that other metal species (usually not 
observed in late B-type stars) are also present in the atmo- 
spheres of these hot BHB stars, but based on our data we 
cannot draw any conclusions about this. Very weak lines 
were used only for the star having the best spectrum. Our 
results thus do not change if we remove the abundance 
determinations from those weak lines. 

In general, our abundance measurements in NGC 1904 
give very similar results to those obtained previously by 
Behr et al. (1999) in M13, which, on the other hand, has 
a very similar metallicity to that of NGC 1904. The very 
small scatter in Fe abundances among the warmer stars 
is particularly intriguing. While it could be the result of 
small statistics, it might indicate that the stars warmer 
than 11000 K observed in NGC 1904 all rotate at very 
similar rates. In fact, Michaud (2002, private comunica- 
tion) estimates that small differences of about 2-3 km/s 
in rotational velocity could generate sizeable discrepan- 
cies in the radiative levitation and gravitational settling 
processes. Unfortunately, the accuracy of our rotational 
measurements in NGC 1904, with errors of the order of 
4-5 km/s, do not allow us to explore this possibility. 
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